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^ . Recently, the question of a relevance of quantum chaos has been discussed in 

applications to quantum computation |^. Indeed, according to the general 
^T) [ approach to closed systems of finite number of interacting Fermi-particles (see, 

Q ■ e.g. 1^, |]), as the interaction between quibits increases a kind of chaos is 

expected to emerge in the energy spectra and structure of many-body states. 
^ . Specifically, the fluctuations of energy levels and components of the eigenstates 

Q-i. turn out to be very large and they are described by the random matrix theory, 

g I Clearly, if this happens in a quantum computer, it may lead to a destruction 

■ of the coherence (due to an internal decoherence inside many-body states) 

required for quantum computations. It is important to stress that the quantum 
chaos occurs not only in the systems with random interactions, but also for 
purely dynamical interactions. In the latter case, the mechanism of chaos is 
^ I the non-linear two-body interaction represented in the basis of non-interacting 

c3 [ particles. 

Numerical analysis |I| of the simplest model of a quantum computer (2D 
model of 1/2-spins with a random interqubit interaction J) shows that as the 
number, L, of qubits increases, the chaos threshold Jcr decreases as Jcr oc 1/L. 
Consequently, it was claimed that the onset of quantum chaos is a real danger 
for the quantum computers with large L ^ 1. On the other hand, in ^ is was 
argued that in order to treat this problem properly, one needs to distinguish 
between the chaotic properties of stationary states and the dynamical process 
of quantum computation. 

Below, we report our theoretical and numerical results for a realistic model 
of quantum computer, described in |[. We consider both stationary and 
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dynamical approaches to the model in the region of a non-selective excita- 
tion which prepares a homogeneous superposition of = 2^ states needed to 
implement both the Shor and the Grover algorithms. 

The model describes a 1-dimensional chain of L interacting 1/2-spins in the 
constant magnetic field B^, subjected to a sum of p = 1, P time-dependent 
rectangular pulses of a circular polarized magnetic field rotating in the x, y- 
plane |]^, |[. Each of the pulses has the amplitude frequency Up, phase ipp, 
and is non-zero during the time Tp = tp^i — tp. The Hamiltonian has the form, 

A;=0 n>k 



1 ^ L-1 / 

^ p=l k=0 \ 



-iupt-iipp J- _|_ ^iupt+iipp j+ \ Q\ 



where the "pulse function" 6p(t) is 1 during the p-th pulse. The quantities 
Jk,n are the constants of Ising interactions between two qubits, Uk are the 
frequencies of the spin precession in the -B"^— field, and flp is the Rabi frequency 
corresponding to the p-th pulse. The operators = ± il^ are defined by 
the relations 1^'^'^ = (l/2)cr^'^'^, the latter being the Pauli matrices. 

The Hamiltonian for a single pulse can be written in the coordinate system, 
rotating around z-axes with the frequency Up. Thus, for one pulse the model 
is described by the stationary Hamiltonian (below, ipp = TT/2,Qp = Q,h'p = 
u). We mainly study the nearest-neighbor interaction {N -interaction) between 
qubits for the dynamical case, Jk,n = J 5n,k+i^ and when all Jk,k+i are random. 
In contrast to the model with homogeneous magnetic field [|I|, we consider a 
constant gradient magnetic field with linear dependence on the position of the 
A;-th qubit, 6k = \uJk+i — ^^k\ ^k = ak, with Qp ^ Jk^n ^ ^ uJk- Thus, 
for the dynamical A^— interaction the Hamiltonian reads, 

h=j:[- dJl + nil] - 2^ E h = ^k- y. (2) 

fc=0 A:=0 

For this Hamiltonian we have developed a theory |Q ^ which predicts two 
transitions which depend on the interaction J. The first transition was called 
in the delocalization border which corresponds to the transition to weak 
chaos for, 

J > J,, ^ — . (3) 

By weak chaos we mean a kind of randomness in many-body states, together 
with the absence of the Wigner-Dyson (WD) distribution P{s) for the spacings 



2 



between energy levels of the Hamiltonian (^. The latter distribution is a strong 
evidence of quantum chaos in the energy spectra of chaotic quantum systems. 
It typically emerges above the delocalization border [^]. Instead, the form of 
P(s) in our model is very close to Poisson, which is known to occur in integrable 
systems. Our analytic approach allows one to explain this unexpected result 
by showing that, indeed, the model (@) is close to the integrable one, even in 
the case of a completely random A^-interaction [Q. 

The estimate turns out to be very different from that obtained in [jl]] 
for a homogeneous magnetic field. Indeed, according to (0), the (weak) chaos 
border is independent of the number of qubits. Therefore, a magnetic field with 
a constant gradient strongly reduces the unwanted effects of quantum chaos. 
Numerical data show that one needs to have a relatively weak interaction, 
J <^ J(,r, in order to avoid large errors in the structure of many-body states, 
which appear as a result of weak chaos. 

Another unexpected analytical prediction which is confirmed by the nu- 
merical data, is that the delocalization border (^) remains the same for the 
case when all qubits interact with each other with random interactions, Jk^n- 
However, in this case, the delocalization border coincides with the on- 
set of strong chaos. The latter is characterized by strong (almost Gaussian) 
fluctuations of the components of eigenstates, and by a WD-distribution for 
P{s). Theoretical analysis shows that the mechanism of this transition to 
strong chaos is related to a strong overlap of energy bands in the spectra of 
the Hamiltonian (0). 

We have also studied the errors that arise when preparing the uniform 
many-body state from the ground state. For this, we computed the evolution 
of the wave function in the model (|I|), during one pulse with = 7i/2. Without 
the interaction, J = 0, and at the absence of the magnetic field gradient, at the 
end of the pulse all components of the wave function are the same, ijj^ = 1/\/N. 
The interaction causes in some errors which can be characterized by the am- 
plitude, T] = {\\ipn\ —i^nDn, and the phase, = (arctan(/m'?/'„/i?e-?/'„))„, where 
(...)„ means the average over different n components. Numerical data show 
that the errors decrease with an increase of Q as rj oc Q'"^, and (p in 
agreement with simple analytical estimates. As one can see, the delocalization 
border does not influence the errors. This means that weak chaos is not im- 
portant for this kind of evolution ( "non-selective" ) of our system. Indeed, this 
evolution lasts only a short time (r = tt/Q) compared with the inverse distance 
between nearest levels inside the energy band. Therefore, when the bands are 
non-overlapped, weak chaos does not influence the dynamics. On the other 
hand, as fl decreases, the bands start to overlap which strongly increases errors 
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in the wave function. 
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